10 PRINT "EXAMPLE 7.5-2 JUMP IN A TRAPEZOIDAL CHANNEL": PRINT
20 DEFINT I: DEF FNC (DY) = Q ^ 2 * (B + 2! * M * DY) - G * (DY * (B + M * DY)) ^ 3'EQ.(12.5.7)
30 DEF FNFM (DY) = .5 * B * DY ^ 2 + M * DY ^ 3 / 3! + Q ^ 2 / (G * DY * (B + M * DY)) 'EQ.(12.4.3)
40 READ G, Q, B, M, GAM: DATA 9.806,16.,4.,.4,9806.
50 YMAX = 16!: YMIN = 0!: PRINT : PRINT "G,Q,B,M,GAM="; G; Q; B; M; GAM
60 FOR I = 1 TO 15:  YC = .5 * (YMAX + YMIN)
70 IF FNC(YC) > 0! THEN YMIN = YC ELSE YMAX = YC
80 ' PRINT YMAX; YMIN; YC
90 NEXT I: PRINT "Y1,YC="; .5 * YC; YC
100 Y1 = .5 * YC: YMIN = YC: YMAX = 3! * YC: FM = FNFM(Y1)
110 FOR I = 1 TO 15: Y2 = .5 * (YMAX + YMIN)
120 IF FNFM(Y2) - FM > 0! THEN YMAX = Y2 ELSE YMIN = Y2
130 ' PRINT "YMAX, YMIN, Y2="; YMAX; YMIN; Y2
140 NEXT I
145 A1 = Y1 * (B + M * Y1): A2 = Y2 * (B + M * Y2): V1 = Q / A1: V2 = Q / A2
150 LOSS = (V1 ^ 2 - V2 ^ 2) / (2! * G) + Y1 - Y2: POWER = GAM * Q * LOSS / 1000!
160 PRINT "Y1,Y2, VI, V2, LOSS, POWER=";
170 PRINT USING "###.### "; Y1; Y2; V1; V2; LOSS; POWER